DOE/ER/40427-6-N96 
ADP-96-24/T223 



Determination of flavor asymmetry for E ± 
by the Drell-Yan process 

"vo . 

Os . Mary Alberg a ' b , Ernest M. Henley b ' c , Xiangdong Ji d * and A.W. Thomas 6 

ON 



Oh' 
CO 



X 



a Department of Physics, Seattle University, Seattle, WA 98122, USA 
Department of Physics, Box 351560, University of Washington, Seattle, WA 98195-1560, USA 



Massachusetts Institute of Technology, Cambridge, MA 02139, USA 



c National Institute for Nuclear Theory, Box 351550, 
University of Washington, Seattle, WA 98195-1550, USA 

>. 

Center for Theoretical Physics, Laboratory for Nuclear Science and Department of Physics, 

ON 
O 

\Q . e Department of Physics and Mathematical Physics, 

and Institute for Theoretical Physics, University of Adelaide, Adelaide, S.A. 5005, Australia 

Flavor asymmetries for the valence and sea quarks of the S 1 * 1 can be ob- 
tained from Drell-Yan experiments using charged hyperon beams on proton 
and deuteron targets. A large, measurable difference in sea quark asymme- 
tries is predicted between SU(3) and pseudoscalar meson models. The latter 
predict that in S + , u/d < 1/2, whereas the former predict u/d ~ 4/3. Esti- 
mates of valence quark asymmetries based on quark models also show large 
deviations from SU(3) predictions, which should be measurable. 
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Parton distributions contain a wealth of information concerning the non-perturbative 
structure of hadrons. While most attention has so far been paid to the distributions of 
the nucleon, through the Drell-Yan process one can also get information on the parton 
distributions of the hyperons. We shall be concerned with the predictions of the distibutions 
for the £ hyperons in various models, including: (1) a quark and meson model, (2) an SU(3) 
model and (3) a model involving the coupling of octets of mesons and baryons (without 
regard to mass). In particular, a comparison of £ and proton structure functions can 
be used to differentiate between these models. We show that measurements of the Drell- 
Yan process for ^ scattering on protons and neutrons (in deuteron targets) allow one 
to extract information on the parton distributions of the £ and that the expectations 
are quantitatively quite different for (1), (2) and (3), above. In particular, we find that 
measurements of the u/d ratio in the £ + are very sensitive to the model chosen, and that 
a quark-diquark model of the hyperon predicts large SU(3) violations in the valence quark 
distributions. 

One of the surprises in the structure of the proton is that the sea appears to have a flavor 
asymmetry, an excess of d compared to u Although the experimental results could 

imply some violation of isospin, this appears to be less likely, and we interpret them as an 
SU(2)q flavor asymmetry in the sea @]. Thus, the d excess in the proton is expected to be 
reflected in an excess of u in the neutron; isospin symmetry would be broken if this were 
not the case. The evidence for flavor asymmetry in the proton sea is based on analyses of 
deep inelastic muon scattering |],[| and Drell-Yan processes 0]. One explanation that has 
been offered is that the excess of d over u is due to the Pauli exclusion principle A 
more likely explanation, in our view, is that offered by Thomas and colleagues Henley 
and Miller |]13[ , and others p4|-pl[| , namely that the presence of a pion cloud surrounding a 



proton favors d over u because of the excess positive charge of the meson cloud. 

It is interesting to apply these arguments to the strange baryons and to compare them 
with SU(3). Here we focus on the charged £ + (£~), composed of uus (dds) valence quarks. 
Thus the main difference from the p(n) case is the replacement of a valence d(u) quark by an 



s quark. In the following, the quark distribution q(u, d, or s) without subscripts refers to the 
proton, and with a S subscript refers to the S + . The x-dependence of these distributions is 
implied, but not shown explicitly. With neglect of mass effects or under SU(3), 

r = -= = —— = 0.51 i 0.04 i 0.05, (1) 

d s s 

where the experimental ratio || is that obtained for the proton at x ~ 0.18. 

The ratio k = 2s/ (d + u) is a measure of the strange quark content of the nucleon. It 
has been determined experimentally in neutrino-induced charm production |22|-p4| to be in 



the range 0.373^;^ ± 0.018 < k < 0.57 ± 0.09. The CTEQ p| determination of parton 



distributions from global QCD analyses of experimental data uses the value k = 0.5, from 
which 



^ 1 -£ = 0.25, (2) 



so with d ~ 2u from Eq. 1, 



u + d 2 



-1 = ^ = 1^0.75. (3) 



We also find 



^ = I « 0.38 (4) 

St. d 

from the same analysis. 

For the model of quarks surrounded by light pseudoscalar mesons, the S + has an ex- 
cess of d over u (and the opposite for the £~) in contradistinction to the SU(3) predic- 
tion. If we neglect higher masses, then the T^{uus) will have components A°(uds)ir + (ud), 
S° (uds)n + {ud) , T, + (uus)tt (^[dd — uu}), or p(uud)K°(ds); similarly a T,~(dds) can be 
A°(uds)n~(du), Y, c \uds)Ti~ (du) , £~ (dds)7r c \^\dd — uu]), or n(udd)K~(us). Thus there 
is a clear enhancement of d for S + and u for E~ and we expect fs < 0.5. 

We also consider an SU(3) model in which a baryon is composed of octets of baryons 
and mesons. We use the SU(3) isoscalar factors and representation matrices given by the 
Particle Data Group f26|. For 8i — ► 8 (g) 8 



N 



9i 



[3A% -Nrj- 3EK - AK] 



(5) 



and for 8 2 — ► 8 <g> 8 



N 



92 



VSNn + V3Nr] + - V3AK 



(6) 



The standard D and F couplings are related to g\ and by D = ^Wgi and F = ^gi, so 



V 



V8D 



pUl + a)n° + V3 (a - r? j - v^l + ct)mr+ + \/2(a - 1)S + ^ C 



+ (1 - a)T?K + - V3 ( a + - ) AK+ 



(7) 



f . ^uu — dd r- ( 1\ uu + dd — 2ss 
p\(l + a) — + \/3 f a - ~ 



v/2 



v/6 



2(1 + a)nud 



+v / 2(a - + (1 - a)Y?us - V3 [a + - j Aus 



with a = F/D, and which leads to relative probabilities, averaged over x, 



d 



2 r 

9 
2 

9 
4 



9ol + Qa + 1 



9a 2 + 18a + 13 



s ~ - 



18a 2 - 12a 



Then 



r = 



m w s 1 + 6F/D + 9{F/Df 



d s s 13 + 18F/D + 9(F/D) 2 



and 



dy 



r s = 



8-12F/D + 18{F/Df 



u + d m e + s e 7+ 12F/D + 9(F/D) 2 ' 
With a = 0.6, consistent with a recent analysis [p^| , we obtain for the proton 



r = .29, f, = .42 



(8) 



(9) 
(10) 

(11) 



(12) 



(13) 



(14) 



which differ significantly from the experimental result (f = .51) and parameter (f s = n/2 = 
.25). We also show in Table I the prediction of this model, fs = 0.54, which, like the meson 
cloud model, disagrees with the SU(3) expectation of 4/3. 
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Deviations from SU(3) predictions are also expected for the valence quark distributions 
in E + On the basis of SU(3) symmetry we expect 

r E = — « - ss 0.57(1 -x). (15) 
me u 

The functional form is taken from a fit by CDHS [J28| , and agrees with the latest parton 
distribution analysis of CTEQ |25j within 20%, which is adequate for our calculations. 

We find that for x > 0.2, quark models predict valence quark flavor asymmetries in 
the S + that are greater than the SU(3) result, e.g. by a factor of 3.4 at x — 0.7. Our 
approach to estimating valence quark distributions in the S + is based on a quark-diquark 
model initiated at Adelaide |7|,|29| which has led to the study of charge symmetry violation 



in the nucleon ||30|| . It was found that the dominant contribution to the structure function 
q(x) in the valence region comes from a state in which the two spectator quarks are in their 
ground states. The effective mass of this diquark state will deviate from 3/4 of the nucleon 
mass (in the MIT bag model, 2/3 in the constituent quark model) because of the hyperfine 
interaction. The mass difference between the two spin states of the diquark leads to spin 
and flavor dependence of q(x) pi|| . Let x q represent the most probable momentum fraction 
carried by the quark q, and x qq represent the most probable momentum fraction carried by 
the diquark. Then the peak in q(x) can be estimated from 

Xq ~\- Xqq 1, Xq 1 Xqq ~ 1 , (16) 

TUB 

in which m qq and m# are the diquark and baryon masses, respectively. For the nucleon, 
the N — A splitting leads to m qq = 650 MeV in the spin singlet state, and m qq = 850 MeV 
in the spin triplet. Then in the proton, d(x) peaks at Xd ~ 0.10, whereas u(x) peaks at 
x u ~ 0.31 - at the scale appropriate to the model. After QCD evolution, these estimates 
are in reasonable agreement with recent parton distribution analyses p5| . 

These same arguments may be applied to the S + . In this case the diquark uu must be 
in a spin triplet, so from the A — £ splitting, m uu S3 850 MeV, and ss(^) peaks at x s ~ 0.28. 
This is close to the value found for u(x) in the proton, so we set ss(x) ~ u(x)/2 (the factor 
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of 2 comes from normalization). To estimate the «s distribution we note that the su diquark 
mass is increased by 180 MeV because of m S) and with the hyperfine splitting, m su 900 
MeV in the singlet, leading to a peak x u (S — 0) ~ 0.24 - i.e., a "harder" distribution, like 
that of u(x) in the proton - and m su m 1050 MeV in the triplet, leading to x u (S — 1) « 0.10 
- a "softer" distribution like in the proton. Since the singlet and triplet diquark states 
are equally probable, we approximate uj;(x) ~ d(x) + u(x)/2. Then 

rs ~ m s ~ 2(rf + w/2) ~ 1 + 2d/u ~ 1 + 1.14(1 -x)' [ ' 

In the valence quark region, this ratio is considerably in excess of that predicted by SU(3), 
as can be seen in Fig. 1, in which we plot the ratio R 

(rs) OT( 3) ~ [1 + 1.14(1 -x)]0.57(l-x) 1 ) 

There are a number of ways that these arguments and models can be tested for S 
hyperons. The most practical to us appears to be in terms of the Drell-Yan cross sections 
for T^p and S ± n (i.e. d) - e.g., in the inclusive reactions T^p — > l + l~X, where P 1 are muons 
or electrons and X is unmeasured. Beams of appear to be adequate for this purpose, 
but n contamination will lead to problems which need to be overcome. 

We first consider the determination of sea quark flavor asymmetry for E ± . We find 
that extracting the ratio f^(x) = us(x)/d^(x) for the S + depends on the known ratios 
r(x) = u{x)/d{x) and f(x) = u/d in the proton. The former is well-determined, and the 
recent determination of the latter has been discussed above. Ratios involving s in the E ± 
cannot be tested easily because they involve second order annihilations (ss on s), terms 
which we neglect because the present accuracy of Drell-Yan measurements is insufficient to 
be sensitive to them. 

Drell-Yan cross-sections are proportional to the products q(x)q(x'), weighted by the 
product of the quark charges, and summed over contributions from beam and target. We 
neglect sea-quark - sea-quark collisions, which would contribute below the likely level of 
accuracy of the experiment. We assume isospin reflection (charge) symmetry: u(x) = d n (x), 



u{x) = d n (x), = d^-(x), uz+(x) = dz-(x), and se+(x) = s^-(x). In the following 

equations, q(x) represents valence quarks and q(x) represents sea quarks. The valence quark 
normalizations are: / u(x) dx = 2 and / d(x) dx = 1. 

Consider the Drell-Yan process for EiV. Let <r(£iV) represent the cross-section for inclu- 
sive dilepton production 

<Pa(EN -*l + l~X) 8na 2 rjr . N _ 9r . . . , rv ^ 

(T(EiV) = s 1 d ^- dy J - = l^K(x^ x N ) J2 e?{?< Wft(^) + P ^ iV] (19) 

with M the mass of the dilepton pair and y/r = Mj \fs. The factor K(xs, %n) accounts for 
higher-order QCD corrections. If the cm. rapidity y ~ 0, then x^ ~ x^ ~ x, and 

87TQ; 2 4 1 

cx(E + p) w — -=iT(a;){-[w(x)w s (x) + w s (x)w(x)] + -\d{x)d Ti {x) + s E (x)s(x)]}. (20) 
9y t y y 

Then by charge symmetry 

87TCt 2 1 _ 4 

cr(E n) w — -^.fr(x){-[w(x)w s (x) + u E (x)w(x) + s E (x)s(x)j + -d(x)d T ,(x)} . (21) 
yy r y y 

We also find 

8ira 2 4 - 1 - 

o-(£ + n) — -=K(x){-[d(x)uz(x) + Us(x)d(x)\ + ~[u(x)dv(x) + s E (x)s(x)]}, (22) 

yyr y y 

and again by charge symmetry 

8im 2 1 - 4 - 

<r(E p) ~ — -=Ar(x){-[d(x)M S (x) + u E (a;)d(a;) + s s (x)s(x)] + -u(x)cfe(x)}. (23) 
yyr y y 

As we note below, if K(x) is known, and all four cross sections are measured, u%, ds, u-s, 
and can be determined. The uncertainties in K(x) can be factored out by taking ratios 
of cross sections; two independent ratios can be constructed. We first define a ratio R'(x) 
determined from the Drell-Yan cross-sections so as to eliminate all unknowns except for 
f s (x) 

= [a(£+p) - q(S-n)] + f(x)[a(£-p) - a(E+n] 
l J " [a(E+p)-a(E+n)]+A[cr(E- p )-a(E-n)} ' 1 J 

and use Eq. ^ - Eq. ^ to write R'(x) in terms of the ratios fs(x), r(x) and f(x): 
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_ r s (x)[r(x) - r(x)] - [1 - r(x)r(x)] 
K {X) ~ 5[r(x) - 1] • { ' 

Thus for r(x) ~ 2 and f(a;) ~ 0.5, R'(x) ~ 0.3fs(x). 

If is known, c?e(x) can be determined directly from the cross sections: 

j M 27y^ [g(S+p) - gg+rQ] + 4[a(S-p) - agTn)] 
S{ ) A^a*K{x) [u(x) - d(x)] ' 1 ' 

and se(^) can De determined from the cross sections and s(x): 

27y^ [g(S+n) - 4a(S-p)] - r (.t) [g(S+p) - 4o-(S~n)] 
EW 87ra 2 K(x) s(x)[r(x)-l] ' 1 ' 

(Recall that, because of the higher mass of the strange quark, we expect s%(x) to peak at a 

larger x than d(x) -c.f., Eq. 11.) 

Quark models with a meson cloud predict the sea quark distributions q(x); they also 

predict that the difference D = x[d(x) — u{x)\ peaks at x ~ 0.1 p~4|-f2"0[|. On the basis of 



meson cloud models^, the distributions of sea quarks in the £ may differ somewhat from 
those in the nucleon due to the presence of kaons; this may shift the maximum of D to 
somewhat smaller values of x. Nevertheless, the region < x < 0.2 should be a good one in 
which to determine fs- 

We believe that the measurement of R' should be possible to within w 20% and this is 
sufficient to establish the preponderance of (J over u in the E + , as predicted by the octet and 
meson cloud models. From Eq. ^5], an error, e, in the measurement of R' leads to an error 
of approximately 3e in f^. If fs were found to be < 0.5, together with the known value of 
f pa 0.51, this measurement would help to reinforce the necessity to include pseudoscalar 
mesons in quark models of baryons. 

To measure valence quark asymmetries we consider the Drell-Yan process for S + and S~ 
on isoscalar targets - with cross sections cr+ and <r_, respectively. We fix xs to be above 0.3 
so that valence quarks in the hyperons dominate. Then from Eq. |TP|- Eq. ^ with x = x-% 
and x' = xtv, 



l We are undertaking a calculation of the S 1 * 1 sea quark distributions. 
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and 



Rirry 2 A A 2 

a + = a(Z+A) » + + -[a E (s) (28) 



a = *{1TA) « ^^(x,x'){^ S (x)[«(x') + <iV)] + ^(^(x')]}- (29) 



We approximate if (x, a;') by an average K, and integrate over x' in the nucleons, so that 

J dx'a + (x',x) = *^^K{^u s (x)[u + d) + ^[s s (x)s)}, (30) 
and similarly for cr_, with q = J dx' q(x'). Then 

~ f dx' a+(x' , x) ~ 4ux(x)(d + u) + 2sj:(x)s ~ 4 + /tr E ' 

We again use the CTEQ f25[ value, k = 0.5, and evaluate R v for both the SU(3) prediction 

for r E (Eq. |15D and for our quark model (Eq. [L7]). In Fig. 2 we plot the ratio D 

R v (quark model) 
D[X) " fl.(SU(3)) ' (32) 

We note that the predicted asymmetry exceeds the SU(3) prediction by about 10% at x = 
0.5, increasing to 20% at x = 0.75. An accuracy of ~ 5% should be possible for these 
integrated cross sections. Thus these measurements will test the SU(3) violations in the 
valence quark distributions of E predicted by quark models. 

In summary, there are substantial differences expected between the valence and sea par- 
ton distributions associated with several models of hyperon structure. We have seen that 
Drell-Yan experiments based on existing hyperon beams should be capable of testing these 
ideas. The substantial violations of SU(3) flavor symmetry in the valence distributions are 
probably the easiest to test as they require only an isoscalar target and a semi-integrated 
cross-section. However, the enormous interest in the underlying cause of the flavor asym- 
metry of the proton sea should also make the tests of sea quark distributions an important 
priority as well. 
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TABLES 



TABLE I. Sea quark asymmetries 



meson cloud 



SU(3) 



octets 



experiment 



r = — 

— d 



ry, = 



u+d 

dT, 



theory ref 



< 0.5 



4/3 



0.29 
0.42 
0.54 



0.51 
= 0.25 
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FIGURES 

FIG. 1. Plot of the ratio R, Eq. 18, as a function of x. 
quark model, Eq. 17, and as found from SU(3), Eq. 15. 



We also show rs, as obtained from the 



FIG. 2. Plot of D, Eq. 32, as a function of x. 
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